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ABSTRACT 

We present a cross-correlation analysis to constrain the faint galaxy contribu- 
tion to the cosmic X-ray background (XRB). Cross-correlating faint optical galaxy 
catalogues with unidentified X-ray sources from 3 deep ROSAT fields we find that 
B < 23 galaxies account for 20 ± 7% of all X-ray sources to a flux limit of 
S(0.5 — 2.0keV) = 4 x 10 _15 ergs _1 cm -2 . To probe deeper, galaxies are then cross- 
correlated with the remaining unresolved X-ray images. A highly significant signal is 
obtained on each field. Allowing for the effect of the ROSAT point-spread function, 
and deconvolving the effect of galaxy clustering, we find that faint B < 23 galaxies di- 
rectly account for 23±3% of the unresolved XRB at IkeV. Using the optical magnitude 
of faint galaxies as probes of their redshift distribution, we find evidence for strong 
evolution in their X-ray luminosity, parameterised with the form L x cx (1 + z) 3,2±1 °. 
Extrapolation to z = 2 will account for 40 ± 10% of the total XRB at IkeV. The na- 
ture of the emitting mechanism in these galaxies remains unclear, but we argue that 
obscured and/or low luminosity AGN provide the most plausible explanation. 

Key words: galaxies: evolution - galaxies: active - X-rays: general - X-rays: galaxies 
diffuse radiation 



1 INTRODUCTION 

The nature and origin of the cosmic X-ray background 
remains a major unsolved problem. The most significant 
progress has been made in the soft X-ray band below 3keV 
since the launch of high resolution imaging satellites such as 
Einstein and ROSAT. By resolving as many sources as pos- 
sible in the deepest ROSAT exposures (Hasinger et al. 1993) 
up to 70% of the 0.5 — 2 keV XRB has now been resolved 
into discrete sources. 

Using a survey of 7 deep (21-57ks) ROSAT fields we 
have detected over 400 X-ray sources above a 4a thresh- 
old to an approximate flux limit of 5(0.5 — 2.0 keV) ~ 
4x 10 -15 ergs -1 cm -2 . By optical spectroscopy we then iden- 
tify the optical counterparts to these X-ray sources. This 
technique has shown that broad-line QSOs directly account 
for at least 30% of the total 0.5 - 2keV XRB flux (Shanks 
et al. 1991). As a larger sample of QSOs became estab- 
lished, detailed studies of the QSO X-ray luminosity function 



(Boyle et al. 1994) and the source number count distribu- 
tion (Georgantopoulos et al. 1996) have shown that QSOs 
are unlikely to contribute more than 50% of the XRB at 
IkeV. The existence of a new X-ray emitting population 
was postulated. 

A further problem in explaining the XRB with QSOs 
is the shape of their X-ray spectra. QSOs show relatively 
steep X-ray spectra with indices of F = 2.2 ±0.1 while the 
1 — 10 keV XRB has a significantly flatter spectrum with 
T = 1.4 (Gendreau et al 1995). Deep ROSAT surveys have 
begun to resolve a population of harder X-ray sources at 
the faintest flux limits (Hasinger et al. 1993, Vikhlinin et 
al. 1995, Almaini et al. 1996). There have been suggestions 
that this hardening may be due to a change in the intrinsic 
properties of QSOs at high redshift or the effect of interven- 
ing photoelectric absorption (Morisawa & Takahara 1990, 
Vikhlinin et al. 1995) but recent work by Almaini et al. 
(1996) found no significant change in the 0.5 — 2keV X-ray 
spectra of QSOs with flux or redshift. The spectral harden- 
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ing is due to another, largely unidentified X-ray population. 
Many of these unidentified X-ray sources appear to be asso- 
ciated with faint galaxies, with implied X-ray luminosities 
typically 100 times higher than normal field galaxies. In par- 
ticular, a number of individually identified narrow-emission 
line galaxies show spectra significantly harder than QSOs, 
more consistent with the spectrum of the XRB (Almaini et 
al. 1996). Similar results were obtained by Carballo et al. 
(1995) and Romero- Colmenero (1996). 

The nature of the X-ray emitting mechanism in these 
galaxies remains unclear. This will be discussed further in 
Section 4, but it should be emphasized that many of these 
objects could contain low luminosity or obscured AGN (see 
eg. Comastri et al. 1995) . In this sense the entire XRB could 
still be due to 'AGN' rather than two distinct classes of X- 
ray source. Whatever the source of X-ray emission, it has 
been known for some time that objects classified as galaxies 
could make a significant contribution to the XRB. Using 
Einstein observations of local galaxies, Giacconi & Zamorani 

(1987) extrapolated the observed X-ray to optical flux ratios 
to the faintest (B < 27.5) optical counts of Tyson et al. 

(1988) and found that normal galaxies could produce at least 
13% of the XRB at 2keV. Griffiths and Padovani (1990) 
used the observed evolution in the 60^im luminosity of IRAS 
galaxies and the local X-ray to infrared ratios to estimate 
that IRAS and starburst galaxies seen to high redshift could 
produce 10-30% of the 0.5-3keV XRB. 

Other studies have probed the galaxy contribution using 
the two-point cross-correlation between the hard (2-10 keV) 
XRB and various optical and infra-red galaxy catalogues 
(Lahav et al. 1993, Miyaji et al. 1994, Carrera et al. 1995, 
Refregier et al. 1997). These galaxy catalogues were gener- 
ally fairly shallow (z < 0.1) but by correcting for the effects 
of clustering and extrapolating the measured volume emis- 
sivity to z~5 it was shown that 10-30% of the XRB could 
be explained by sources associated with galaxies. However 
because of the unknown error in such an extrapolation and 
the uncertain evolutionary properties of galaxies it is clearly 
desirable to probe more distant objects directly. The clear- 
est evidence that fainter galaxies could be important con- 
tributors to the XRB came from the study of Roche et al 
(1995). Using the improved sensitivity and angular resolu- 
tion of the ROSAT satellite they were able to probe fainter 
limits and much smaller angular scales (~ 15"). By cross- 
correlating faint galaxy catalogues with X-ray sources they 
obtained a 2 — 3u detection implying that B < 21 galax- 
ies account for ~ 5% of the X-ray sources. However they 
obtained a more significant (~ 5a) signal by removing the 
X-ray sources and cross-correlating with 3 deep (21-49ks) 
images of the unresolved XRB. The results implied that 
galaxies to a limit of B = 23 directly contribute ~ 17% 
of the IkeV XRB. A simple extrapolation to B < 28 showed 
that galaxies could account for ~ 30% of the XRB or pos- 
sibly more with evolution. Similar results have since been 
obtained by Roche et al. (1996b) and Soltan et al. (1997). 

In this paper we perform an independent test of these 
results on 2 new deep (~50ks) ROSAT exposures and for 
the first time attempt to measure the evolution in the X-ray 
emissivity of faint galaxies with redshift. The deepest field 
(GSGP4) from the Roche et al. (1995) work is included, but 
this time the effect of clustering is analysed more rigorously. 
First we will perform a cross-correlation of the unidentified 



Field 


RA 


DEC 


N H 


Exposure 


BJS855 


10 h 46 m 24 s 


-00°21'00" 


2.9 ±0.4 


57147s 


BJS864 


13 h 43 m 4 3S 


-00°15'00" 


2.6 ±0.3 


52466s 


GSGP4 


00 h 57 m 29 s 


-27°38'13" 


1.8 ±0.3 


48995s 



Table 1. Summary of deep ROSAT fields, with coordinates in 
J2000 and galactic column density in 10 20 atom cm -2 . 

X-ray sources with faint galaxies. We will then probe deeper 
into the remaining unresolved XRB. To deconvolve the ef- 
fect of galaxy clustering we will make use of the method 
developed by Treyer & Lahav (1996) (hereafter TL96), ex- 
tending their formalism to take account of the ROSAT point 
spread function. We will then consider the cross-correlation 
signal as a function of magnitude to obtain an estimate of 
the evolution of X-ray emissivity with redshift. We assume 
the cosmology q — 0.5 and A = throughout. 

2 CROSS-CORRELATING X-RAY SOURCES 
WITH FAINT GALAXIES 

2.1 Data and method 

In this paper we use 3 deep ROSAT PSPC exposures with 
optical identifications from the X-ray source catalogue of 
Shanks et al. (1997). We refer to this catalogue paper for 
full details of the X-ray source detection and optical follow- 
up observations. To a flux limit of 5*(0.5 — 2.0 keV) ~4x 
10 _15 ergs _1 cm -2 , approximately 50% of the X-ray sources 
were identified as QSOs by optical spectroscopy (Boyle et al. 
1994). Many of the remaining X-ray sources had no obvious 
QSO candidate within their error box, but appeared to be 
associated with faint galaxies on deep photographic plates 
(Georgantopoulos et al. 1996). Due to the high surface den- 
sity of galaxies at faint magnitudes however (~10000 deg -2 
at B<23, Metcalfe et al. 1995) and the relatively large X- 
ray error circle (~25" FWHM) many of these will be chance 
associations. If we look deep enough we can find a galaxy 
counterpart to any X-ray source. We will therefore adopt 
a statistical approach to determine the fraction of galaxies 
among the X-ray sources. 

The 3 deepest ROSAT fields from the survey will be 
used in this analysis, including the GSGP4 field analysed by 
Roche et al. (1995). A summary of these exposures is given 
in Table 1. The optical galaxy catalogues were obtained from 
COSMOS plate scans of Anglo- Australian Telescope photo- 
graphic plates, calibrated with magnitudes and a star-galaxy 
classification flag. This optical data was previously used for 
the galaxy clustering and number count analyses of Steven- 
son et al (1985) and Jones et al. (1991). Further details of 
the plate reduction and zero point magnitude calibration are 
given in Jones et al (1991). Astrometric solutions to convert 
the plate co-ordinates into R.A. and Dec. were already in 
existence from the UVX QSO survey of Boyle et al. (1990) 
and Boyle, Jones & Shanks (1991). A second transforma- 
tion was then required to convert R.A. and Dec. into the 
(x,y) ROSAT pixel coordinates. This was performed using 
approximately 10-15 bright stars and UVX QSOs with defi- 
nite counterparts on each ROSAT image. These gave 6 coef- 
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Figure 1. Showing the positions of 8714 galaxies with magni- 
tudes B < 23 on the BJS855 field after transformation to X-ray 
coordinates (scaled to degrees). 




ficient astrometric transforms with rms positional errors of 
between 7.1 and 8.4 arcseconds, as expected given the an- 
gular resolution of the ROSAT PSPC. On Figure 1 we show 
the resulting positions of B < 23 galaxies on one of the 
fields (BJS855) in X-ray coordinates. The plate edges and 
larger plate holes are clearly visible, where bad data has 
been flagged around bright objects or in some case where 
plate defects were detected in the plate measuring process 
(Jones et al. 1991). 

In calculating the cross-correlation of optical positions 
with X-ray sources the number of galaxies were counted in 
successive annuli of 5" width around each X-ray source. This 
distribution was then compared with the counts obtained by 
placing 50,000 random points over each of the field areas, 
taking care to avoid the plate edges and "holes". A two 
point cross- correlation function C xo {9) of X-ray and optical 
sources can then be defined as: 



N xo (6i)N r 



(1) 



N xr {0i)N„ 

where, N r is the total number of random points, N 
the number of optical plate sources being correlated and 
N xo (8i) and N xr (8i) give the number of X-ray/optical and 
X-ray/random pairs respectively. 

2.2 Results 

Removing only the X-ray sources identified with QSOs or 
galactic stars we cross-correlate the positions of the remain- 
ing 149 unidentified 4a sources with the objects classified 
as galaxies on photographic plates. These optical galaxies 
were split into a bright sample with B < 21 and a fainter 
sample with 21 < B < 23. The results are shown in Table 2, 
where we display the number of galaxies found within 30" 
of an X-ray source compared to the distribution expected 
by chance. A 30" radius approximately corresponds to the 
3a ROSAT PSPC positional error. 

Considering first the cross correlation with brighter B < 
21 galaxies, on each field we find an excess compared to a 
random distribution. Overall we find 52 source-galaxy pairs 
compared to 31.5 expected by chance, amounting to a 3.7cr 




150 200 



Figure 2. (a) The cross-correlation function C xo (8) of unidenti- 
fied X-ray sources from 3 deep fields with B < 21 mag galaxies. 
On (b) we repeat the analysis with fainter 21 < B < 23 mag 
galaxies. Note the change in scale compared to (a). Error bars 
are due to y/~N statistics. 



rejection of the null hypothesis that the two distributions 
are not correlated. This would indicate an excess of ~ 20.5 ± 
7.2 B < 21 galaxies around X-ray sources. On Figure 2(a) 
we display the source-galaxy cross-correlation in 15" radial 
bins. 

To probe deeper, we repeat the cross-correlation with 
the ~5000 fainter galaxies in the magnitude range 21 < 
B < 23. We find 326 galaxies within 30" of an X-ray source 
compared to 278 expected. This represents a 2.9<r rejection 
of the hypothesis of no correlation and indicates an excess 
of ~ 48 ± 17 of these fainter galaxies around X-ray sources. 
The cross-correlation function is displayed in Figure 2(b) for 
the 3 fields combined. 

We must now make a correction for the effect of galaxy 
clustering. Since galaxies are clustered on the scales probed 
by this analysis this can enhance the excess of source-galaxy 
pairs. We therefore apply a simple correction based on 
LO gg (8), the angular autocorrelation function of faint galaxies 
(from Roche et al 1996a). Defining AN' gx (8) to be the ob- 
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IV/T a crn 1 1~ 1 1 ri o T"cincrr> 


Field 


Pairs; <r SO" 


n, V r*t t> rl 

J_J A |J l_ L I <^ LI 


rT \~srr"f*w tM crn i fir*n n cc* 




GSGP4 


22 


13.2 


2.4cr 


B < 21 


BJS855 


19 


12.3 


1.9a 




BJS864 


11 


5.93 


2.1a 




Total 


52 


31.5 


3.7a 




GSGP4 


123 


104.9 


1.8a 


21 < B < 23 


BJS855 


101 


85.0 


1.7a 




BJS864 


102 


88.1 


1.5a 




Total 


326 


278 


2.9a 



Table 2. Showing the number of faint galaxies found within 30" of unidentified X-ray sources on 3 deep fields compared to the number 
expected by chance. The galaxies are split into those with B < 21 and fainter galaxies with 21 < B < 23. 



served excess of source-galaxy pairs measured within a given 
angle 6, this will be enhanced relative to the true number of 
X-ray emitting galaxies AN gx (6) such that: 



AN gx (0) = 



AN' gx (0) 



1 + G; o )AN gg (8) 



(2) 



where AN gg (6) is the excess of galaxies within angle 6 
compared to a random distribution, obtained by integration 
of u> gg {9), and Gtot is the total number of galaxies in the 
entire image. For the brighter B < 21 galaxies this gives a 
25% clustering enhancement within a radius of 30". At 21 < 
B < 23 the effect is slightly more significant, accounting for 
31% of the excess signal. 

Thus overall, allowing for the effects of clustering, we 
have evidence that ~ 48 ± 17 of the 149 unidentified X-ray 
sources are due to galaxies with B < 23, strengthening the 
initial findings of Roche et al (1995). Hence from a total of 
236 sources on 3 deep fields we have shown that faint (B < 
23) galaxies contribute ~ 20 ± 7% of the total source counts 
to a flux limit of S(0.5 - 2.0keV) ~4x lO^ergs^cm- 2 . 

However when we scale this 20% galaxy contribution 
by the median flux of the unidentified X-ray sources (~ 7 x 
10 _15 ergs _1 cm -2 ) we obtain a total contribution of 1.2 ± 
0.4 x 10~ 9 erg s _1 cm ~ 2 sr _1 which accounts for only ~ 4% 
of the total XRB. We will therefore attempt to probe much 
deeper in Section 3 by cross-correlating with the remaining 
unresolved XRB. 



ages used in the cross-correlation are extracted from only 
the 0.5 — 2keV data. Data from periods of high particle 
background were also removed from the analysis, excluding 
approximately 10% of the photons when the Master Veto 
Rate was above 170 counts s _1 (Plucinsky et al. 1993). Ad- 
ditional data was available for the BJS855 and BJS864 fields 
from serendipitous pointings obtained from the ROSAT 
data archive, offset 6 and 7 arcminutes respectively from 
the original field centres and providing an additional ~ 20ks 
of data. This data was processed separately before a mo- 
saic of the final images was produced. The 4a sources were 
then removed from the data using a range of exclusion radii 
suitable for excluding 98% of the photons, depending on the 
off-axis angle of the source. This radius varied from 30.6" 
for a source on axis to 75.6" for a source at a radius of 20 
arcminutes. Inspection of the resulting images revealed that 
only one source showed evidence for being extended and 
could not be successfully removed in this way. This source, 
GSGP4X:032, was identified as a distant (z = 0.56) galaxy 
cluster. A 70% larger radius was required to successfully 
remove this source. Finally, only the central 16 arcminutc 
regions were used in the cross-correlation analysis since the 
sensitivity of the PSPC drops off rapidly beyond 20 arcmin- 
utes (Hasinger et al. 1992). For the BJS855 and BJS864 
fields, where deeper exposures have been obtained by adding 
additional archive data, only regions of the final image lying 
within 16 arcminutes of both field centres are used. 



3 CROSS-CORRELATING THE UNRESOLVED 
X-RAY BACKGROUND WITH FAINT 
GALAXIES 

3.1 The data 

In this section we extend the cross-correlation technique to 
probe the remaining unresolved component of the cosmic 
X-ray background and thus investigate the origin of X-ray 
emission beyond the limit of significant source detection. 
By performing the analysis with different magnitude slices 
of the galaxy sample we may also be able to deduce the 
redshift evolution of X-ray emissivity from faint galaxies. 

The X-ray images were obtained from deep (~ 50ks) 
ROSAT observations reduced using the STARLINK Asterix 
X-ray data reducing package. Since the X-ray background 
below 0.5 keV is dominated by galactic emission and so- 
lar scattered X-rays (Snowden & Freyberg 1993) the im- 



3.2 Cross-correlation results 

A two point cross- correlation between the unresolved X- 
ray images and the faint galaxy catalogues was obtained by 
counting the number of X-ray photons in successive annuli 
around each galaxy and comparing this with the mean pixel 
intensity. By repeating for every galaxy in the field (avoiding 
holed regions), the cross-correlation is defined by: 



W xg (6i) = 



(3) 



where N x (8i) and N p (6i) are the total number of X-ray pho- 
tons and pixels respectively within annulus i of each galaxy, 
(Nx) is the mean number of photons per pixel and the sum- 
mations occur over all the galaxies in the field. Error esti- 
mates were obtained by splitting each field into 4 quadrants 
and obtaining the cross-correlation for each separately. The 
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18<B<23 




Field 


Photons < 30" 


Expected 


Excess significance 


GSGP4 


8575 


8154 


4.7 


BJS855 


5763 


5490 


3.7 


BJS864 


7771 


7414 


4.1 


Total 


22109 


21058 


7.2 



Table 3. Results of the cross-correlation of 18 < B < 23 cata- 
logue galaxies with the unresolved 0.5 — 2.0 kcV X-ray background 
on 3 deep fields, comparing the number of 0.5 — 2.0 kcV X-ray 
photons found within 30" with a random distribution. 




that faint galaxies are not only significant contributors to 
the resolved X-ray source counts but also to the fainter, as 
yet unresolved component of the X-ray background. In Ta- 
ble 3 we show the number of galaxy-photon pairs within 30" 
of the galaxies compared to the counts expected from a ran- 
dom distribution. The excess X-ray photons close to galax- 
ies suggest > 3.5<t detections on each field assuming Poisson 
statistics. We have clearly confirmed the initial findings of 
Roche et al. (1995). 

In Figure 4 we show W xg for the 3 deep fields com- 
bined, where we have summed galaxy-photon pairs across 
all fields and calculated the expected counts by combining 
the 3 separate ^ N p (6i) (Nx) terms. Errors have been esti- 
mated in the same manner as before by calculating 4 cross- 
correlations from the quadrants of each field and using the 
error on the mean values of W xo - 





18<B<23 




BJS864 









100 200 300 



Figure 3. The cross-correlation function W xg (0) of the unre- 
solved 0.5 - 2kcV X-ray background with faint 18 < B < 23 
galaxies on 3 deep fields. 

error on W xg was then estimated from the error on the mean 
from these quadrants. 

The resulting two point cross-correlations between faint 
18 < B < 23 galaxies and the unresolved XRB are shown 
in Figure 3 in 15" bins. It is clear that a significant positive 
signal has been detected on each field. The errors on the 
first two bins alone (< 30") suggest individual detections 
> 3<t significance. The independent analysis of the GSGP4 
field is seen to be in excellent agreement with the results of 
Roche et al. (1995). Thus we find further strong evidence 



3.3 The effect of galaxy clustering 

We have established a highly significant cross-correlation be- 
tween the unresolved component of the cosmic XRB and 
faint galaxies. Galaxies do not randomly sample the sky 
however and are well known to show evidence of cluster- 
ing and structure on the scales probed by our analysis. It 
is therefore probable that a non-negligible fraction of the 
enhanced signal in the cross-correlation is due to the clus- 
tering of galaxies with each other. In effect, the galaxies will 
correlate with the emission of their neighbours as well. One 
approach (see Roche et al. 1995) is to apply an approximate 
correction using the angular clustering of the observed op- 
tical galaxies. However this will not account for the X-ray 
emission arising from fainter unseen objects (B > 23) which 
may be clustered with the galaxy catalogue. This effect is 
discussed in detail in the work of TL96, where a prescrip- 
tion for modelling these populations is presented. We will 
now apply this formalism to our data, taking care to allow 
for the effect of the ROSAT PSPC point-spread-function 
(PSF). This effect was neglected by TL96. A detailed de- 
scription of the effect of the PSF can be found in Refregier 
et al. (1997). 

3.4 Correlation functions and the effect of the 
PSF 

Under the assumption that the XRB arises from discrete, 
point-like X-ray sources, it has been shown that two terms 
will contribute to the cross-correlation of faint galaxies with 
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100 200 
Separation (arcsecs) 



Figure 4. The total cross-correlation function W xg (8) of the un- 
resolved 0.5 — 2.0 keV X-ray background on 3 deep ROSAT fields 
with 18 < B < 23 galaxies. Also shown is the best fitting model of 
the form given in Equation 14, formed as the sum of the Poisson 
term (dotted line) and a clustering term (dashed line). 



the X-ray background (Lahav 1992, Treyer & Lahav 1996, 
Refregier et al. 1997). 



Vo = Vp + Vc 



(4) 



For the ideal case of a delta-function PSF with cell sizes 
u> — > 0, this function is related to the normalized form of W xg 
(Equation 3) by: 

Vo 



W xg {9) = 



LU 2 IN 



(5) 



In this expression / is the mean intensity of the unre- 
solved XRB and N gives the number of catalogue galaxies 
per steradian. The r\ v term comes from the direct contribu- 
tion from the catalogue galaxies themselves. This is known 
as the Poisson term, which in the idealized case only exists 
at zero lag (6 — 0). For cells of solid angle lo this term is 
directly related to the X-ray intensity AI g contributed by 
the galaxies such that: 



Vp 



Pa( z ) j i n 
47rrF 



(6) 



where we define p g (z) as the observed volume emissivity 
of the galaxy population, dv(z) is the volume element per 
unit angle and ri is the luminosity distance. 

The second term n c arises from the clustering of the X- 
ray sources with the galaxy population and is known as the 
clustering term. Treyer & Lahav (1996) have calculated the 
theoretical angular cross-correlation of a galaxy population 
with the X-ray background. They find that the clustering 
term ij c obeys a simple power law of the form: 



where the amplitude A xg can in principle be used to 
evaluate the X-ray volume emissivity due to galaxies given 
a specific model of their redshift distribution, clustering and 
luminosity evolution (see Section 3.5). TL96 applied this 
simple form to the measurements of Roche et al. (1995) to 
obtain an estimate of the galaxy contribution to the XRB. 
As we will demonstrate below, fitting this functional form 
directly to measurements of W xg (9) isn't strictly valid due 
to the effect of ROSAT PSF. This has the effect of smearing 
out the Poisson term to make a significant contribution at 
non-zero lag. We will discuss how to deconvolve the Poisson 
and clustering terms in Section 3.6. First we will use the 
formalism of TL96 and a specific model for the galaxy pop- 
ulation to relate the amplitude of the clustering term A xg 
to the X-ray volume emissivity. 



3.5 Modelling the galaxy population 

We assume that any diffuse component of the 0.5 — 2.0 keV 
XRB is negligible in comparison with the source component 
and model the observed volume emissivity of the source pop- 
ulation using an evolutionary parameter q such that: 



Pxrb(z) = p s (z) = po(l + z) q 



(8) 



Thus the intensity of the XRB per unit solid angle is 
given by: 



I — U) 



Ps(z) 
4nrf(z) 



dv(z) 



(9) 



We will further assume the spatial cross-correlation of 
galaxies with the X-ray sources £(r, z) to be the same as the 
auto-correlation of faint galaxies with themselves. 

Next we must characterise the properties of the galaxy 
population at a given redshift. Deep spectroscopic sur- 
veys (eg. Glazebrook et al 1995) have measured the galaxy 
N(m, z) to optical magnitudes of B ~ 24. We hope to model 
even deeper than this however. Fainter galaxies lack red- 
shifts, so we will use the latest evolution model of Roche et 
al. (1996a). This model is consistent with current spectro- 
scopic results (eg. Cowie et al. 1996) and reproduces the faint 
number count distribution. In general this model predicts a 
higher redshift distribution for faint blue galaxies than the 
parametric form of Efstathiou (1995) used in TL96. 

Next we must model the clustering properties of the 
galaxy population and its evolution with redshift. We will 
adopt the standard form (Peebles 1980): 



ar,z) = (l + z)-^ (l-X 



(10) 



vc(0) = A xg e 



1-7 



(7) 



where the proper coordinate r is the spatial separation 
between the sources and e models the clustering evolution. 
We will adopt a stable clustering model (e = 0) and a cor- 
relation radius of r = 4.4 /i _1 Mpc (Loveday et al. 1995). 
Using these parameters and the galaxy evolution model de- 
scribed above, Roche et al. (1996a) were able to fit the ob- 
served galaxy angular clustering to a magnitude of B = 27. 

To interpret our observations we will use the framework 
developed by TL96, who have calculated the theoretical an- 
gular cross-correlation of a galaxy population with the X-ray 
background. The detailed calculations are somewhat cum- 
bersome and will not be repeated here. However they find 
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that the clustering term ij c (6) obeys a simple power law 
(Equation 7) where the amplitude is given by: 



47T 

The function H 1 is given by: 



dx(l + x 



-7/2 



(11) 



(12) 



Assuming the standard 7 = 1.8 leads to #(7) = 3.68. 
Defining a global evolution parameter K = 7 — e + q ~ 5, the 
function /(N) takes the form: 



f 

J B, 



N(m, z)dm. 



(13) 



where B m in and B max represent the magnitude range 
of the galaxy catalogue. 

Thus we may in principle re-arrange Equation 11 to 
obtain a value for the local X-ray volume emissivity p us- 
ing our chosen description of the galaxy population and our 
observed X-ray parameters I and A xg . As discussed above 
however, measuring A X g is non-trivial since we must decon- 
volve the effect of the Poisson term. This will be discussed 
in the next section. 

In order to estimate the evolutionary parameter H we 
need a value for the parameter q which describes the evolu- 
tion in the observed X-ray emissivity (Equation 8). We will 
attempt to measure this evolution directly in Section 3.7. 

3.6 Deconvolving the Poisson and clustering 
terms in W xg (9) 

We define the PSF of the ROSAT PSPC to be the function 
ip(6) normalised to unity over the whole sky. We adopt the 
parametric form given by Hasinger et al. (1995). Our images 
are extracted from the central 16 arcminute region of the 
PSPC and we will use the mean weighted PSF over this 
radius, assuming a mean photon energy of IkeV. 

The effect of the ROSAT PSF is to smear out the X- 
ray emission over neighbouring pixels, and hence the Poisson 
and clustering terms in the cross-correlation, derived assum- 
ing a delta function PSF, will be smeared out compared to 
Equation 5 to produce the measured W xg (6) as follows: 



W xg (8 12 ) = 



IN2-KUJ 2 



d0i 



d^i 



d^ 2 



/ dfisV(M7?(6>le 

J s 



(14) 



The df2i and d^l 2 integrals are evaluated over cells C\ 
and C2 separated by #12. A rigorous derivation of such equa- 
tions can be found in Refregier et al. (1997). The final in- 
tegral is evaluated over the whole sky, and is effectively a 
2-dimensional convolution of the point spread function with 
the idealized correlation function ??(#). 

Although the two terms 7/ c and r/p both contribute to the 
observed W xg (9) they are not independent. By Equation 6 
the Poisson term r\ v can be evaluated from the X-ray volume 
emissivity, which in turn is related to the clustering term rj c 
via Equations 8 and 11. Therefore r]{6) in the above integral 
may be replaced by: 




1 1.5 
Redshift 



Figure 5. Galaxy redshift distributions predicted using the PLE 
model of Roche et al. (1996a) for B < 23 galaxies split into a 
'bright' and 'faint' sample. 



V(0) 



Air 



/(i± 

J z 4tt 



4nrf 



-dv 



(15) 



We must therefore evaluate this functional form of 
W xg (0) numerically and then fit to our observations to ob- 
tain an estimate of the amplitude A xg . In order to evaluate 
this function, however, we must assume a value for the evo- 
lutionary parameter q. We will now attempt to measure this 
evolution directly. 



3.7 Constraining evolutionary parameters 

The largest uncertainty in previous estimates of the galaxy 
contribution to the XRB has been the assumed redshift evo- 
lution (eg. see Lahav et al. 1993, Roche et al. 1995). Treyer 
& Lahav (1996) proceed under the assumption that the lu- 
minosity ratio Lx / L opt remains constant at all redshifts. To 
study this assumption further, we consider cross-correlating 
with successively fainter magnitude slices of the galaxy pop- 
ulation. Since fainter galaxies, on average, will probe higher 
redshifts it should in principle be possible to constrain the 
redshift evolution of the galaxy population. Modifying Equa- 
tions 11 and 13 we can obtain theoretical cross-correlations 
of the unresolved XRB with galaxies between the magnitude 
limits [ml, m2]: 



dAx9 ~ MAN J z dZ{1 



+ 2) r 



N 1-7 



N(m,z)dm (16) 



We can normalize this by the observed cross-correlation 
from the full (18 < B < 23) galaxy sample, giving: 



A N I' f™ 2 

SA * 3 = JcMan J dz{1 + z)Hr1 "' J iV(m ' 2)dm (17) 

Thus by performing the cross-correlation with galaxies 
in different magnitude ranges we may be able to constrain 
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Figure 6. The cross-correlation function W xg (0) of galaxies with 
the unresolved 0.5 - 2.0 keV XRB from three ROSAT fields. The 
galaxies arc split into a bright and faint sample. 



the evolutionary parameter K and hence the evolution in 
X-ray emissivity. 

The galaxy catalogue was therefore split into "bright" 
(18 < B < 21) and "faint" (21 < B < 23) subsets. The 
redshift distributions for these samples are shown in Fig- 
ure 5, as predicted from the models of Roche et al. (1996a). 
Two separate cross-correlations with the unresolved XRB 
are then carried out. The results are displayed in Figure 6. 
The brighter sample clearly show a higher overall correlation 
amplitude. To compare with a range of theoretical predic- 
tions (in order to constrain the evolution) we evaluate the 
integral in Equation 17 to obtain expected amplitude A xg 
as a function of the evolutionary parameter K. These are 
displayed in Figure 7. Clearly we expect the relative cross- 
correlation amplitudes to change dramatically depending on 
the redshift evolution. With stronger evolution, the ampli- 
tude for the brighter (more local) sample will decrease. 

To constrain K we therefore evaluate 5A xg for both 
galaxy samples. We fit models of the form given by Equa- 
tions 14 and 15. A complication arises since the functional 
form we fit requires a value for q in advance (to evaluate 
the Poisson contribution), so we start with the assumption 
that q = and proceed iteratively. In practice the relative 
amplitudes of the Poisson and clustering terms are fairly in- 
sensitive to the assumed value of q and vary by only ~ 10% 
with q in the range q £ [0, 3]. Nevertheless we start with a 
value q = 0, use the measured amplitude A xg to constrain 
the evolutionary parameter ft (and hence q) and then re- 
peat the process. Within a couple of iterations we converge 
on values for the amplitudes A xg for both magnitude slices. 

With 8 measured in radians we obtain amplitudes of 
5. 15 ±0.35 x 10~ 5 for the bright galaxies and 2.8 ±0.2 x 10~ 5 
for the fainter sample. These amplitudes were obtained by 
fitting to W xg (0) for 6 < 2 arcminutes and the errors are 
derived from the variance in the amplitudes from 3 inde- 
pendent fields. The resulting estimates for the evolutionary 



Figure 7. Showing the expected relationship between the ampli- 
tude of the cross-correlation function W xg and the evolutionary 
parameter K for brighter 18 < B < 21 galaxies (solid upper line) 
and the fainter 21 < B < 23 galaxies (lower dashed line). The 
shaded region displays the la error bounds on the measured am- 
plitude for the brighter galaxies and the corresponding range of 
K. 



parameter K are shown in Figure 7. The shaded area shows 
the la confidence region on the observed amplitude for the 
brighter galaxies and the corresponding preferred values of 
N. For the fainter sample however the weaker dependence of 
A xg on K does not allow us to provide any useful constraints 
on the evolutionary parameter. This is expected however 
since, by construction, Equation 17 is normalized by the 
amplitude of the full 18 < B < 23 dataset and therefore 
dominated by these fainter galaxies. 

The errors on 8A xg suggest values of K in the range 
— 1.37 < K < 0.58. This global parameter was defined 
K = 7 — e ± g — 5 where q describes the evolution of the 
X-ray emissivity as a function of redshift (Equation 8) and 
7 and e characterise the clustering properties of the galaxy 
population (see Equation 10). We have assumed that the 
latter two quantities are relatively well defined in compari- 
son with the parameter q and adopt the values 7 = 1.8 and 
e = 0. This value of e allows the galaxy evolution model of 
Roche et al (1996a) to simultaneously fit the faint galaxy 
number counts and the galaxy clustering. This leads to val- 
ues of q = 2.82 ± 0.98 and thus a best estimate for the 
evolution of the X-ray emissivity due to galaxies: 



p x (z) = p a {\ ± z) 



(18) 



Using a spectral index of a x ~ 0.6 (Almaini et al. 1996) 
we can obtain our best estimate for the X-ray luminosity 
evolution of faint galaxies: 



(1 + *) 



3.22±0.98 



(19) 



which is very similar to the X-ray evolution of AGN. 
Similar forms have been obtained by considering the lu- 
minosity function of individually identified X-ray luminous 
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galaxies at much brighter fluxes (Griffiths et al. 1995, Boyle 
et al. 1995a). 

Since the N(m,z) of Roche et al. (1996a) places the 
faint blue galaxy population at relatively high redshift, it 
might be argued that the strong evolution we have detected 
here is due to our particular choice of galaxy evolution 
model. We have therefore repeated the analysis using the 
model of Efstathiou (1995), following TL96, which assumes 
a much more local population of faint galaxies. Although the 
mean redshifts are lower, this model uses a much smaller cor- 
relation radius (r = 2 h~ Mpc) with comoving clustering 
evolution (e = —1.2) rather than the stable model we have 
used in this paper (e = 0). The end result is a very sim- 
ilar galaxy-galaxy angular autocorrelation. Even with this 
model we still obtain evidence for strong evolution in X-ray 
luminosity (L x oc (1 + z) 3 - 5±1 - 2 ). 



~i — 1 — 1 — 1 — 1 — 

A R] ue (B-R<1.5) 
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3.8 The X-ray volume emissivity of faint galaxies 

With an estimate for the X-ray evolution we may now esti- 
mate the volume emissivity due to faint galaxies using Equa- 
tion 11. First we require the amplitude of the observed angu- 
lar cross-correlation function W xg (9) (Figure 4). As before, 
fitting the functional form given by Equations 14 and 15, we 
obtain (measuring 6 in radians): 

A xg = 3.60 ± 0.27 x 10~ 5 (20) 

The mean intensity of the unresolved XRB from the 3 
deep ROSAT fields used here is / = 1.59 ± 0.2 x 10 -8 erg 
s _1 cm ~ 2 sr _1 and N = 2.64 x 10 7 galaxies sr _1 in the 
magnitude range 18 < B < 23. Putting all of this together 
we may now obtain the local X-ray volume emissivity via 
Equation 11: 

p ~ 3.02 ± 0.23 x 10 38 /iergs _1 Mpc -3 (21) 

It should be emphasized however that this estimate re- 
lies on very specific assumptions about the distribution and 
clustering properties of the galaxy population. In particular 
since p oc r^ 1 (from Equation 11), a less clustered model 
with r < 4.4fa -1 Mpc will increase the required local emis- 
sivity. Such weak clustering would be required in models 
which place most faint galaxies at relatively low redshift (eg. 
Efstathiou 1995). Recent deep spectroscopic work seems to 
confirm the existence of a significant high redshift popula- 
tion however (eg. Cowie et al. 1996). 



3.9 The contribution of faint galaxies to the 
unresolved XRB 

We may now obtain a rough estimate the total contribution 
of the 18 < B < 23 catalogue galaxies to the unresolved 
XRB by integrating the volume emissivity out to the me- 
dian redshift of the sample, z — 0.45. Modifying Equation 
9 gives the following expression for the contribution to the 
sky intensity per unit solid angle: 

= 3.6 ±0.5x10 Vgs 'cm 2 sr 2 



Figure 8. The cross-correlation function W xg (0) of galaxies with 
the unresolved 0.5 — 2.0kcV XRB, separating blue from red galax- 
ies and combining the results from three ROSAT fields. 

This accounts for 23±3% of the unresolved X-ray back- 
ground, where the error is derived entirely from the uncer- 
tainty in the observed X-ray quantities, W xg (9) and I. 

To determine the "total" galaxy contribution to the 
XRB we simply extrapolate Equation 22 to the faintest 
galaxies at arbitrarily high redshifts. This process becomes 
increasingly uncertain as we extend the redshift distribution 
beyond the limit of our sample at B = 23, but nevertheless 
by integrating to z = 2 we formally obtain ~ 80 ± 20% of 
the unresolved 0.5 - 2.0 keV XRB, which is ~ 40 ± 10% of 
the total XRB intensity. Given the uncertainties already in- 
herent in this procedure it will probably suffice to say that 
galaxies can produce a significant fraction of the XRB at 
least as high as the contribution from QSOs. 

3.10 Red and blue galaxies 

To constrain the colour and type of galaxies producing the 
cross-correlation signal we use R band plates to separate the 
sample into red and blue subsets, dividing at B — R = 1.5. 
The cross-correlations were then carried out separately for 
each dataset. The results (Figure 8) show no significant dif- 
ference in the cross-correlation amplitudes, suggesting that 
galaxies of all morphological types are contributing to the 
signal. Blue galaxies at these magnitudes are known to be 
more weakly clustered however, with autocorrelation ampli- 
tudes lower by a factor of ~ 4 at these magnitudes (Roche et 
al. 1996a) and hence the enhancement in W xg {9) due to clus- 
tering will be significantly lower for the blue subset. Thus 
the typical X-ray luminosities may be higher for the blue 
galaxies (by ~ 50%), but a more careful analysis is required 
to prove this since the N(m, z) distribution for blue and red 
galaxies may be significantly different. 
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4 FUTURE PROSPECTS 

Having established the existence of a population of highly 
luminous X-ray galaxies, the next intriguing question is the 
origin and nature of this activity. One explanation could be 
the presence of a large, hitherto undetected population of 
low luminosity Seyfert 1 galaxies. This would certainly ex- 
plain the cross-correlation signal. So far however the small 
number of individually identified X-ray luminous galaxies 
have shown only narrow optical emission lines. High reso- 
lution optical spectroscopy suggests a mixture of starburst 
and Seyfert 2 activity (Boyle et al 1995b, McHardy et al. 
1997) but in many cases the classification was ambiguous. 
Iwasawa et al. (1997) argue that the ambiguous line ratio di- 
agnostics may be due to the presence of an obscured active 
nucleus surrounded by significant starforming activity in the 
host galaxy, with the AGN producing the hard X-ray emis- 
sion. Spatially resolved optical spectroscopy would be able 
to test this hypothesis. Obscured AGN provide a natural 
explanation since the expected X-ray absorption can read- 
ily reproduce the flat spectra of the XRB (Comastri et al. 
1995, Madau et al. 1994). At least two of the brighter X-ray 
galaxies identified in Almaini et al. (1996) show clear evi- 
dence of X-ray obscuration. The discovery of a high redshift 
counterpart to this population adds further credence to this 
possibility (Almaini et al. 1995). Explaining the XRB with 
an AGN population may be difficult however since QSOs are 
strongly clustered and may violate the upper limits of the 
isotropy of the XRB (Georgantopoulos et al. 1993, Soltan & 
Hasinger 1994). Little is known about the clustering char- 
acteristics of low luminosity and obscured AGN however. 

Starburst activity alone is unlikely to explain the X-ray 
emission in most of these galaxies. There have been sug- 
gestions that massive X-ray binaries formed in the wake of 
star formation in early, low metallicity epochs may provide 
a strong source of hard X-ray flux (Griffiths and Padovani 
1990, Treyer et al. 1992). Such models have difficulty in ex- 
plaining the XRB however, since ASCA observations of star- 
burst galaxies reveal very soft X-ray spectra which cannot 
contribute significantly to the hard XRB (eg. Della-Ceca et 
al. 1996). Many individually identified X-ray galaxies are 
also several orders of magnitude brighter than any known 
starburst galaxy. A third explanation has emerged recently 
based on advection dominated accretion onto quasars at 
low accretion efficiency (Di Matteo & Fabian 1997), but the 
enormous black hole masses required present significant dif- 
ficulties for this model. 

We conclude that the obscured AGN hypothesis pro- 
vides the most plausible explanation for the origin of the X- 
ray background and the bulk of the X-ray activity in these 
galaxies. Infra-red spectroscopy may provide a conclusive 
test of this model by allowing us to see through the ob- 
scuring dust to detect broad emission lines in moderately 
obscured AGN. If the nuclei are very heavily obscured then 
the infra-red emission may also be obliterated, but in this 
scenario it is difficult to produce the soft X-ray flux we ob- 
serve with ROSAT without a very large ~ 10% scattered 
component or some additional source of soft X-ray photons. 
Spectrapolarimetry may then allow us to detect broad emis- 
sion lines scattered around the obscuring medium and into 
our line of sight. Further observations are clearly required 



before we can claim to understand the nature of this X-ray 
activity. 



5 SUMMARY AND CONCLUSIONS 

By cross-correlating faint galaxy catalogues with unidenti- 
fied X-ray sources a strong (4.2<r) signal was detected in- 
dicating that individual galaxies with magnitude B < 23 
account for 20 ± 7% of all X-ray sources to a limiting flux 
of ~ 4 x 10 _15 ergs _1 cm -2 in the 0.5 - 2.0 keV band. This 
builds on the results of Roche et al. (1995) who found a 
significant signal in cross-correlation with brighter B < 21 
galaxies, attributing ~ 6% of the X-ray sources to these 
brighter objects. Scaling the 20% galaxy fraction by the me- 
dian flux of the unidentified X-ray sources leads to ~ 4% of 
the total XRB intensity. 

To probe deeper we cross-correlate with individual pho- 
tons in the remaining unresolved XRB images. Significant 
signals were obtained on all 3 deep ROSAT images, each 
of similar amplitude, independently confirming the results 
obtained by Roche et al. (1995). To translate these cross- 
correlations into a total fraction of the unresolved XRB 
a specific description of the galaxy population is adopted, 
modelling the evolution, number density and clustering 
properties of the faint blue galaxy population using the for- 
malism developed by TL96. We modify this formalism to 
allow for the effect of finite cell sizes and the effect of the 
ROSAT point-spread function. By comparing the theoret- 
ical XRB cross-correlation with the observed signal an es- 
timate for the local X-ray volume emissivity was obtained 
at p ~ 3.02 ± 0.23 x 10 38 /iergs _1 Mpc~ 3 . Extrapolation to 
high redshift (z = 2) suggests that faint galaxies can account 
for ~ 40 ± 10% of the total XRB at IkeV. These estimates 
have a strong dependence on the assumed distribution and 
clustering properties of the faint galaxy population. 

When the optical galaxy catalogue is separated into 
blue and red subsets, dividing at B — R = 1.5, no signif- 
icant difference is found in the cross-correlation with the 
XRB. This would suggest that a mixture of galaxy colours 
and morphologies contribute to the observed signal. Given 
that red galaxies are more strongly clustered at these mag- 
nitudes, this may also indicate that bluer galaxies are in- 
trinsically more X-ray luminous. 

To constrain the evolution of X-ray emissivity with red- 
shift, separate cross-correlations were carried out with two 
magnitude slices of the galaxy population. The resulting 
difference in amplitude suggests that faint galaxies evolve 
strongly with redshift such that: 

L x oc (I + £) 3 ' 22±0 - 98 (23) 

which represents the first evidence that the X-ray emis- 
sion from faint blue galaxies evolves as strongly as AGN. 
Similar results have been obtained by analysing the brightest 
narrow emission-line galaxies emerging from deep ROSAT 
exposures (Griffiths et al. 1996, Boyle et al. 1995a). 

Combined with recent findings which suggest that X-ray 
luminous galaxies have hard X-ray spectra, it now seems es- 
tablished beyond doubt that a population of faint galaxies, 
or some processes associated with them, are emitting vast 
amounts of X-ray radiation which may finally solve the puz- 
zling origin of the X-ray background. Obscured AGN models 
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provide a plausible explanation, but the true nature of these 
X-ray luminous galaxies remains unclear. 
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